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Abstract 

I present a spatial analysis of the galaxy distribution around the cluster CI 0024+17. The basic 
aim is to find the scales where galaxies present a significant deviation from an inhomogeneous 
Poisson statistical process. Using the generalization of the Ripley, Besag, and the pair correlation 
functions for non-stationary point patterns, I estimate these transition scales for a set of 1,000 
Monte Carlo realizations of the CI 0024+17 field, corrected for completeness up to the outskirts. 
The results point out the presence of at least two physical scales in this field at 31.4" and 1 12.9". 
The second one is statistically consistent with the dark matter ring radius (~ 75") previously 
identified by Jee et al. (2007). However, morphology and anisotropy tests point out that a clump 
at ~ 120" NW from the cluster center could be the responsible for the second transition scale. 
These results do not indicate the existence of a galaxy counterpart of the dark matter ring, but 
the methodology developed to study the galaxy field as a spatial point pattern provides a good 
statistical evaluation of the physical scales around the cluster. I briefly discuss the usefulness of 
this approach to probe features in galaxy distribution and N-body dark matter simulation data. 
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1. Introduction 

Galaxy clusters are the largest virialized structures in the universe. These systems are com- 
posed of three components behaving differently during collisions: galaxies, hot gas and dark 
matter. While the hot intergalactic gas presents strong eletromagnetic interactions during clusters 
enconteurs, both dark matter and galaxies are predicted to be collisioness, see e.g. 11811 . Hence, 
clusters are important laboratories for studying the dynamics of the different kinds of matter in 
the universe. Indeed, galaxy clusters are themselves the result of several mergers of smaller 
groups of galaxies, see e.g. HlOfl . a process still ongoing in many systems like, for instance, the 
so-called 'Bullet Cluster' (IE 0657-56) [EH], and other clusters like CI 0152-1357 and MS 



1054, see EMU. Also, cluster CI 0024+17 has been a target of many studies since its discovery 
]. This is an intermediate-redshift system (z=0.395) with both weak and strong lensing 
,171]. There are several indications that this cluster, an apparently relaxe d sy stem, is actually 
the portrait of a collision of two clusters, along our line of sight ||8t |9l |20j |26j ll6|]. The weak 
lensing analysis presented by II 1611 shows a ring-like structure in the projected matter distribution 
at r ~ 75" (~ 0.4 Mpc), surrounding a soft, dense core at r < 50". The authors interpret this 
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substructure as the result of a high-speed line-of-sight collision of two massive clusters ~ 1 -2 Gyr 
ago. An interesting question in this context refers to the behaviour of the galaxy component with 
respect to dark matter. In the case of the Bullet Cluster, CI 0152-1357 and MS 1054-0321, the 
distribution of galaxies follows dark matter quite well \17, (jlHi, 15]. However, the recent work 
of 12111 . based on the 2D distribution of galaxies, suggests that the ringlike structure observed in 
dark matter measurements is not seen in the projected two-dimensional galaxy distribution. 

In the present work, I introduce an alternative approach to probe both in radial and angular 
coordinates the projected galaxy distribution around galaxy clusters, and applied it to the case of 
CI 0024+17. The basic aim is to find the scales where galaxies present a significant deviation 
from an inhomogeneous Poisson statistical process and compare that to the dark matter ring 
radius and the secondary clump present in the field. 



2. Galaxy Sample 

The sample of galaxies around CI 0024+17 used in this work is taken from the catalog of 
Czoske et al. (2001). This corresponds to 295 galaxies in the range 0.37 < z < 0.41, the 
most probable cluster members in the ~ 0.3° x 0.3° field. The center of the cluster was chosen 
to be the geometric center of the dark matter ring, see lfl6ll . Distance r for each galaxy was 
calculated from their equatorial coordinates, in arcsecs. Since the completeness of the sample 
varies from >80% at the cluster center to <50% in the outer regions, the catalog was corrected 
through 1,000 Monte Carlo realizations of the sample. This was done following approximately 



the completeness variation map of [8], adopting a procedure similar to [21]. One example of 
such realizations is presented in Figure 1, where black points are taken from the simulation and 
green points are the real galaxies from [8]. The background in this figure is a density map based 
on the void distance from each point with respect to a fine grid ylllj]. This is just introduced to 
illustrate how the simulations fill the empty cells left by the original sample. 



3. Methodology 

3.1. Spatial Point Patterns 

Galaxy distribution around CI 0024+17 is the point pattern we want to probe. In statistical 
theory, a spatial point pattern (SPP) x = {x\, ...,x„] is simply a collection of points x-, in some 
bounded area, and interpreted as a realization of a more generic point process X. Analysis and 
modelling of SPP is a challenging topic in modern spatial spatistics. The mathematical theory 
was first developed to solve various problems where it is sensible to model the locations of events 
as random, lfl3n . Recently, statistical theory has made many developments with respect to SPP 
analysis. One of these advances is related to achievements in the study of nonstationary pro- 
cesses. Statistical stationarity (or homogeneity) means that the distribution of the point process 
is invariant under translations. This implies that the intensity A, i.e. the expected number of 
points per unit area of the SPP, is constant. Nonstationarity, otherwise, allows for deterministic 
local variation of the intensity:/} — > A(u), where u is a given location on the SPP. This is the 
case we are interested in this work, since galaxy distribution in the sampling window around the 
cluster is obviously inhomogeneous. 
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Figure 1 : Simulated + real field on a void distance map. Simulated and real points are in black and green, respectively. 
The background is a density image built from the void distance matrix from each point with respect to a fine grid. 



3.2. Ripley K-function 

In spatial statistics, Ripley K-function lEstl is a classical tool to analyse SPP. It is defined as the 
expected number of events within the distance r of any given point divided by the intensity A. The 
K-function, being a measure of the distribution of the inter-point distances, captures the spatial 
dependence between different regions of a SPP. Mathematically, the inhomogeneous K-function 
is defined as 



— y — 



1{0 < ||k - jc/H < r} | m e X 

) 



(1) 



Thus, A(u)K(r) is the expected total 'weight' of all random points within a distance r of the point 
//, where the weight of a point x, is l/A(xi) yj], and for a Poisson process with intensity function 
A(u), the K-function reduces to K- m j wm — nr 2 . 

A standard computational estimator of K is given by 

Ki„hom(n = — > . > . — ~ ~ - — e(Xi, xy, r), (2) 

area(W) %(Xi)A( Xj ) 

where e(x,, Xf, r) is an edge correction weight, W is the observational window, and A(u) is an 
estimate of the intensity function A(u). 

Although the K-function is the classical tool to probe SPP, modern spatial statistics rarely uses 
K(r), but rather its variant, the L-function as introduced by yj. In the planar case, the L-function 
estimator is defined as 
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The advantage of using L(r) is that it is always proportional to r, and in the Poisson case L(r) - r, 



01311 . Thus, L(r) - r > indicates clustering with respect to a Poisson process. 

In this work, the L-function was applied to the galaxy distribution around CI 0024+17 using 
the estimation algorithms provided in the library spatstat, developed by [3], and running under 
the R statistical package (www.cran.r-project.org). The edge correction used here is implemented 
in the tasks Kinhom and Linhom. It is briefly described in the next. 

3.3. Edge Corrections 

Edge corrections are used to correct biases in the estimation of K(or L)-function due to the 
loss of information near to the boundary of W. In this work, I have applied the following edge 
correction: 

m > r) 

e(Xi,xf,r) - — , (4) 

where bj is the distance from x, to the boundary of the window W. Among other options avail- 
able in spatstat, this weighting function is considered the most computationally and statistically 
efficient when there are large numbers of points, see J2]. 

3.4. Intensity estimator A(u) 

Finally, there remains the question of how to estimate the intensity function A(u). This is 
estimated using a ieave-one-out' kernel smoother, as described in The estimate A(xi) is 
computed by removing x, from the point pattern, then applying a Gaussian kernel smoothing to 
the remaining points, and finally evaluating the smoothed intensity at x, (in this work we set the 
kernel width equal to 30")- This estimator is also included in Kinhom and Linhom. 



4. Analysis of CI 0024+17 

I have used the routine Linhom to probe the galaxy distribution around CI 0024+17. For 
each MC realization of this field, I found the scales where clustering weakens towards an inho- 
mogeneous Poisson process. In Figure 1, regions above the line L(r) = r describes clustered 
distribution of points. To test the null hypothesis of an inhomogeneous Poisson process I also 
plot in this figure the 95% envelopes computed from M independent simulations of this process, 
with the same intensity as the galaxy field, i.e., W\r) for j = 1, ...,M. The upper and lower 
envelopes of these simulated curves are 

L hwer (r) =rrdn£ ( -'V) (5) 

L" pper (r) = max L u, (r). (6) 

If data come from a Poisson process, then L(r) and L (1) (r),...,L (M) (r) are statistically equivalent 
and independent, so the probability that L(r) lies outside the envelopes is equal to 2/(M + 1) 
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by symmetry. This corresponds to the significance level a in the test which rejects the null 
hypothesis. Hence, to compute the 95% envelopes (a = 0.05), we need only 39 MC simulations 
of the Poisson process, see e.g. yj]. 

In Figure 2, we see two different regions where the null hypothesis of a Poisson process is 
rejected at the 5% significance level. The transition scales where the clustered process weakens 
towards a Poisson process are marked with green dashed lines. A single transition around a 
cluster indicates its outer limit, the region where galaxy distribution declines in clustering (or 
intensity) and turns indistinguishble from a Poisson process. Two (or more) transitions indicate 
a recoverage of the clustering process. Any morphological feature could be responsible for that: 
clumps, rings, tails, etc. In the present analysis, about 88% of the MC simulations present two 
scales like those we see in Figure 2. The Shapiro-Wilk's test fails to detect that the distributions of 
1 st and 2 nd transition scales significantly deviate from normal distributions: /Ji =31 .4"; cr = 8.4" 
(p-value=0.3603) and ^ = 112.9"; a = 21.4" (p-value=0.7338). Also, Bartlett's and Kruskal- 
Wallis's tests strongly reject the hypotheses of equal variances and means, respectively, for these 
two distributions (p-value < 2.2 x 10~ 16 in both tests). All these results indicate the presence of 
two distinct physical scales in the CI 0024+17 galaxy field, as we can see in Figure 3. Note that 
the dark matter ring is located at r ~ 75" [16], and separations between this radius and the 1 st and 
2 nd scales are equal to ~ 5.17cn and ~ 1.77o"2, respectively. This suggests that the dark matter 
ring is more likely to be associated to the second region. Also in this figure, it is presented the 
average correlation function for all set of MC simulations (see the small box). Note that the 1 st 
transition scale corresponds just to a narrow diminishment interval in the correlation function, 
while the DM ring radius and the 2 nd transition scale are approximately coincidents with the 
highest peaks of %(r) after ~ 50". Not by chance this scale roughly marks the intersection of the 
two Gaussian distributions. These results indicate that the system has a core at < 50" (which 
agrees with lfl6Tl ) and that there is an extended region after this radius with two significant peaks 
(~ 3cr above £). One of them agrees well with the position of the DM ring radius while the other 
coincides with the mean of the 2 nd transition scale. 

5. Morphology test 

One should be careful not to overinterpret the previous analysis by concluding that the second 
scale we see in Figure 3 implies that the galaxy distribution around CI 0024+17 has a luminous 
counterpart of the dark matter ring. Indeed, this result may have to do with the way observations 
were done. For instance, the completeness map was constructed after a smoothing procedure 
using a Gaussian kernel of 30" [9]. This size is similar to that I have found for the first transition 
scale, but the separation between the two normal distributions is ~2.7 times this smoothing scale, 
so it seems that the width of the kernel is not large enough to invalidate the results. On the 
other hand, the presence of a secondary peak in the distribution about 120" NW of the cluster 
center (see Figure 1 and Figure 4 of |9J]) can be the responsible for the second scale detected 
in data. This can be seen in Figure 4, where I plot the perspective view of galaxy distribution 
also smoothed by a Gaussian kernel of width 30". Letters A and B indicate the first and second 
density peaks in the field. Note that the peaks are very close to each other, although clearly 
distinct. Less distinguishble, however, can be the effect of a spheroidal clump from a ring in 
the spatial analysis I have applied to data. One way to measure such morphologic effects is now 
presented. 

I have built two mock samples made to mimic the cases we want to study. In Model I, we have 
a central cluster with radius equal to /ii = 31.4", plus a clump at ~ 120" NW (radius defined as 
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Figure 2: Besag L function for one realization of the CI 0024+17 galaxy field (solid lines). The Poisson process corre- 
sponds to the line L(r)=r (dashed lines). The 95% lower and upper envelopes (averaged over all the MC realizatios of the 
galaxy field) are presented in dotted lines. Transition scales are marked with vertical green dashed lines 




r (arcsec) 



Figure 3: Normal distributions of the 1 st and 2 transition scales. Vertical dashed lines indicate the position of the DM 
ring radius. Small box: average correlation function for all MC simulations. DM ring, 1 st and 2 nd scales are indicated by 
green dashed lines. 
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Figure 4: Perspective view of galaxy distribution around CI 0024+17 smoothed by a Gaussian kernel of width 30". 
Letters A and B indicate the the first and second density peaks in the field. 
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//i/2) and a Poisson background with same intensity of the CI 0024+17 field. Both the cluster 
and the clump have a Hernquist profile (Hernquist 1990), but the cluster is set to have ~ 4.5 times 
the density of the clump (after the Gaussian smoothing). Model II is the central cluster plus a 
ring at the second scale, 1 12. 9", with width equal to cr 2 = 21.4". The density of the cluster is also 
~ 4.5 times the density of the ring. For each case, 1,000 mock fields were generated. Examples 
are presented in Figure 5, pannels a and c. 

Analysis of these mock fields shows that: (i) in both cases two statistically distinct scales are 
found - Bartlett and Kruskal-Wallis tests reject the null hypotheses of equal variances and means 
at >99% confidence (see resultant normal distributions in Figure 5, pannels b and d); (ii) however, 
while in Model II about 94% of the realizations present two scales, this only happens for 61% of 
the Model I realizations; (iii) the 1 st Gaussian is narrower for Model I {/u\ = 30.1, cr\ = 7.7,//" = 
31.3, erf = 9.7) while the 2 nd Gaussian is narrower for Model II (j/ 2 = \\13,a\ = 29.0,^ = 
120.0, (Tj = 19.3); (iv) however, there is no evidence for difference in the results of Models 1 and 
2 - Bartlett, ANOVA and Tukey tests indicate only a 0.001% probability of difference between 
the samples. 

The morphology test presented here is quite naive in the sense that neither model reproduces 
the complexity of the real galaxy field. Exactly for this reason, it may indicate that the methodol- 
ogy presented in Section 3 is not accurate enough to discriminate different morphologies around 
the galaxy cluster. To improve the results, I have applied a wavelet analysis to the galaxy field 
of CI 0024+17. Here I followed the method of B2411 . Basically, the coordinates of the galaxies 
are convolved by a radial Mexican Hat filter on a grid of N x N pixels. The method leads to a 
multiscale analysis which starts below the biggest scales, where the whole cluster is detected as a 
unique structure, and ends with the smallest scale, when approximately only one galaxy lies into 
the operating area of the wavelet. To obtain the significance level of a given structure I have ap- 
plied the replica procedure, where many simulations are made by drawing independently X, and 
Yi from the the X and Y distributions of the sample studied to destroy small-scale correlations, see 
e-g- lfl9ll . The analysis is then performed on the "random fields" for the same wavelet scale than 
the true field. The significance level of a given structure is estimated from the number of images 
N e xt in which maxima greater than the maxima wavelet coefficients are found among the whole 
set of simulations N se t- Thus, Psl = N e _ tt /N set , see 11911 . In Figure 6 we present the isosurfaces of 
the wavelet coefficients corresponding to the scales a - 64, 32, 16, and, 8, respectively (in pixels 
units for a map of 128 x 128 pixels). The wavelet scale used is a = 64 pixels corresponding to 1 .0 
Mpc. In the top left low resolution image (a = 64 pixels), we see the prominent central region 
of the cluster with some density enhament around it, which could suggest an almost symmetric 
substructure at first sight. However, with increasing resolution, this feature disappears, while at 
scales a — 32 and 64 a core + clump structure emerges. Finally, at scale a = 8 (bottom right 
image) just a central peak remains. For the two intermediate scales I estimate the significance 
level of a model core + clump, having found Psl = 0.023 (a = 32) and Psl = 0.045 (a = 16) 
(for 1,000 MC simulations), which suggests that such a model cannot be rejected at these scales. 

6. Anisotropy test 

Spheroidal clumps and rings (and possibly other shapes) may produce approximately the same 
results after spatial analysis. However, in the case of comparing rings and clumps, there is a 
remaining property to be tested: the anisotropy of the signal around the cluster. This second test 
is now introduced. 
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Figure 6: Morphology test. Pannels (a) and (c) exhibit examples of the distribution of points for the cases cluster+clump 
and cluster+ring, respectively. Pannels (b) and (d) show the result of the spatial analysis: normal distributions related to 
the scales detected for both geometries. 
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For the projected distribution, anisotropy can be probed by the reduced second-order moment 
measure 'K of a point pattern (e.g. Illian et al. 2008 and references therein). In this work, we 
estimate % with using the library spatstat under R statistical package. The routine Kmeasure 
executes the following steps: 

1 . it takes a point pattern, 

2. and forms the list of all pairs of distinct points in the pattern, 

3. then computes the vectors that join the first point to the second point in each pair, 

4. and treats these vectors as a pattern of 'points', 

5. finally applyies a Gaussian kernel smoother to them (keeping the kernel width 30", see 
Sections 4 and 5). 

The algorithm approximates the point pattern and its window by a binary pixel image, in- 
troduces a Gaussian smoothing kernel and uses the Fast Fourier Transform to make a density 
estimate k. The density estimate of k is returned in the form of a real-valued pixel image. The 
% estimator is defined as the the expected number of points lying within a distance r max of a 
typical point, and with displacement vector having an orientation in the range [a,fi\. This can be 
computed by summing the entries over the relevant region, i.e, the sector of the disc of radius x 
centred at the origin with angular range [a,/3]. Hence, we can compute a measure of anisotropy, 
A, as integrals of the form 



Recall that Ripley function (both in first and second order) is used to test the hypothesis that a 
given planar point pattern is a realization of a Poisson process. In the case of vectors instead of 
points, the second-order orientation analysis can be done around the center of the field. A nearly 
flat anisotropy profile would be consistent with a Poisson process A « 0. Thus, the idea here 
is to look for significant anisotropics in integral (7) for angular steps A9 = /? - a. The choice 
of AG is arbitrary, and it was set to be A8 = 15°, based in the following analysis. We take 
1,000 realizations of a controlled sample corresponding to a point pattern given by a Poisson 
distribution plus a Hernquist spheroid (Hernquist 1990) located at ~ 100" and making an angle 
of 45° around the center. To justify our choice of A6 we present in Figure 7 the analyses for the 
case for AO = 5°, 15°, and 30°. For A6 = 5°, we still detect the peak, but now there are secondary 
peaks and a more noisy behaviour for (A). For A9 — 30°, the peak is still there, but now it is 
less significant. This result suggests that in the limit of too small A6 we have a noisy anisotropy 
curve (possibly with false peaks), while in the limit of very large A6 the signal can be completely 
lost. Actually, in the approximate range [12°, 18°] the resultant anisotropy profiles are almost 
indistinguishable, with relative differences < 3%. 

Results of the anisotropy test for the averaged samples are presented in Figure 8. Note that 
Model I (although very simplified) reproduces fairly well the anisotropy profile of real data, 
exhibiting a peak around 110° (or the supplementary angle 290°), a direction corresponding to the 
axis crossing the NW quadrant. On the contrary, Model II has a flatter anisotropy profile. In fact, 
Kolmogorov-Smirnov test is consistent with the null hypothesis that Data and Model I samples 
come from the same distribution. This is not true for the KS test comparing Data and Model 
II samples, and Model I and Model II samples (see p-values in Figure 8). I also have applied 
a Welch two-sample t-test to find if the anisotropy profiles are stochastically higher between 
each other. Again, the results indicate that both Data and Model I have profiles significantly 




(7) 
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Figure 7: Anisotropy profiles of a controlled sample (Poisson + clump at 45°). Solid lines for A8 = 15°; dotted lines for 
A9 = 5°; and dashed lines for A0 = 30°. 

higher than that of Model II. Also, there is not a significant difference between Data and Model 
I profiles (see p-values in Figure 8). This all means that the galaxy field around CI 0024+17 is 
more consistent with an anisotropic cluster+clump distribution instead of a cluster +ring one. 

Morphology and anisotropy tests are not totally conclusive, but they suggest that the second 
scale found in data is more probable to be associated to the clump at ~ 120" NW from the 
cluster center than to a ring of galaxies, although this feature cannot be discarded by spatial 
analysis alone. 

7. Summary and Discussion 

By considering the planar galaxy distribution around CI 0024+ 17 as a SPP, I have applied some 
tools of spatial statistics analysis to probe the existence of an extended structured in this cluster. 
The catalog was corrected through 1,000 MC realizations of original the sample, following the 
completeness map of [8]. Based on the Besag L-function curves I have identified two distinct 
physical scales where the null hypothesis that the data came of a Poisson process is rejected at 
95% confidence. Values found for the two scales present normal distributions: p.\ = 31.4"; cri = 
8.4" and n\ = 112.9";cr 2 = 21.4". Both Bartlett and Kruskal-Wallis statistical tests indicate 
these normals are significantly unequal both in variances and means. This reinforces the idea 
of two density peaks in the CI 0024+17 field. Also, the analysis of the averaged correlation 
function for all MC realizations suggests that at ~ 50" may be the maximum scale for the core 
of the system. This scale is approximately at the intersection of the normal distributions found 
for the two transition scales. The DM ring radius and the 2 nd transition scale are approximately 
coincidents with the highest peaks of £(r) after ~ 50". However, morphology and anisotropy 
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Figure 8: Anisotropy test. The anisotropy profiles for the averaged samples are presented as follows: Data - solid lines; 
Model I - dotted lines; Model II - dashed lines, p-values associated to KS test and Welch t-test are presented in the upper 
left corner table. 



tests based on simulations of 1,000 mock fields for cluster+clump and cluster+ring models are 
not conclusive about interpreting the second scale as a ring of galaxies. However, a wavelet 
analysis of the field points out that a cluster+clump model cannot be rejected at scales a — 32 
pixels (0.50 Mpc) and a — 16 pixels (0.25 Mpc). Finally, the anisotropy test also suggests that 
the second scale found in data is more probable to be associated to the clump at ~ 120" NW 
from the cluster center than to a ring of galaxies. 

Summing up this work, the analyses presented in Sections 4, 5 and 6 are not indicative that 
galaxies follow dark matter around CI 0024+17, but the methodology based on spatial patterns 
examination seems to have good sensitivity to find physical scales in galaxy systems (both in 
radial and angular coordinates). It could turn a useful tool to study galaxy (and maybe other 
astrophysical) systems in general. Spatial analysis (followed by subsidiary morphology and 
anisotropy tests) could be a poweful technique to probe remnant substructures of clusters colli- 
sions using galaxy data alone. Analysing a large dataset of galaxy clusters, we could identify 
candidates to bullet-type clusters by the presence of at least two transition scales. Also, the ap- 
proach can introduce a new estimator for galaxy systems radii based only on the scales detected 
in data [22]. Finally, this method can be directly applied to data from N-body simulations of clus- 
ter collisions. Some sensitivity to initial conditions could be tested not only looking for bumps in 
the density projected along the collision axis, but also with the dark matter particles interaction 
range. This is an interesting test to be done, since a recent study suggests that the ring-like fea- 
ture of CI 0024+17 (and similar cases) would be associated to very improbable initial velocity 
distributions B27I1 . 
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